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ABSTRACT 

Given a large-scale mixture of self-interacting dark matter (SIDM) particles and 
baryon matter distributed in the early Universe, we advance here a two-phase accretion 
scenario for forming supermassive black holes (SMBHs) with masses around ^ lO^M© 
at high redshifts z ( > 6). The first phase is conceived to involve a rapid quasi-spherical 
and quasi-steady Bondi accretion of mainly SIDM particles embedded with baryon 
matter onto seed black holes (BHs) created at redshifts z < 30 by the first generation 
of massive Population III stars; this earlier phase rapidly gives birth to significantly 
enlarged seed BH masses of Mbh,*! 1.4 x 10'^ Mq cro/(l cm^ g-i)(Cs/30 km s"^)^ 
during z ^ 20 — 15, where cto is the cross section per unit mass of SIDM particles 
and Cs is the velocity dispersion in the SIDM halo referred to as an effective "sound 
speed" . The second phase of BH mass growth is envisaged to proceed primarily via 
baryon accretion, eventually leading to SMBH masses of Mbh ^ 10^ M©; such SMBHs 
may form either by z ~ 6 for a sustained accretion at the Eddington limit or later at 
lower z for sub-Eddington mean accretion rates. In between these two phases, there 
is a transitional yet sustained diffusively limited accretion of SIDM particles which 
in an eventual steady state would be much lower than the accretion rates of the two 
main phases. We intend to account for the reported detections of a few SMBHs at 
early epochs, e.g., SDSS 1148-1-5251 and so forth, without necessarily resorting to ei- 
ther super-Eddington baryon accretion or very frequent BH merging processes. Only 
extremely massive dark SIDM halos associated with rare peaks of density fluctuations 
in the early Universe may harbour such early SMBHs or quasars. Observational con- 
sequences are discussed. During the final stage of accumulating a SMBH mass, violent 
feedback in circumnuclear environs of a galactic nucleus leads to the central bulge for- 
mation and gives rise to the familiar empirical Mbh — (^b correlation inferred for nearby 
normal galaxies with df, being the stellar velocity dispersion in the galactic bulge; in 
our scenario, the central SMBH formation precedes that of the galactic bulge. 

Key words: accretion, accretion discs - black hole physics - cosmology: theory - 
dark matter - galaxies: formation ~ quasars: general 



1 INTRODUCTION 

On the basis of various observational diagnostics and nu- 
merous case studies, supermassive black holes (SMBHs) are 
now widely believed to be ubiquitous, particularly at the nu- 



clei of both normal and active galaxies (e.g., Kormendy & 
Richstone 1995; Haehnelt 2004). As the central gravitational 
engines to power most energetic activities of quasi-stellar 
objects (QSOs) or quasars (e.g., Salpeter 1964; Lynden-Bell 
1969; Bardeen 1970), SMBHs dynamically impact the for- 
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mation and evolution of host galaxies (e.g., Silk & Rees 1998; 
Page, Stevens, Mittaz & Carrera 2001; King 2003; Murray, 
Quartaert & Thompson 2005). Their most tantalizing man- 
ifestations are the observed Mbh — Mhuige correlation (e.g., 
Magorrian ct al. 1998; Laor 2001; Haring & Rix 2004) and 
its tighter version — the Mbh oc a"^'^ correlation for both 
active and normal galaxies (e.g., Gebhardt et al. 2000; Fer- 
rarese & Merritt 2000; Tremaine et al. 2002), where Mbh 
is the black hole mass, Mbuige is the galactic central bulge 
mass and uj, is the stellar velocity dispersion in the galactic 
bulge. 

Given substantial progresses in probing and under- 
standing the basic physics of SMBHs as well as galaxy forma- 
tion and evolution over past several decades (e.g., Lynden- 
Bell 1969), much still remain to be learned and explored 
because of the somewhat speculative nature of the subject 
to a certain extent. During the extensive Sloan Digital Sky 
Survey (SDSS), the newly reported SMBH with a mass of 
Mbh ~ 3 X 10^ Mq in the quasar SDSS 1148+5251 (Fan 
et al. 2003; Willott et al. 2003) at a high redshift z = 6.43 
particularly highlights the outstanding mystery of the rapid 
BH mass growth in the early Universe and reveals inconsis- 
tency with the local Mbh — cri relation (e.g., Walter et al. 
2004).^ 

For a sustained Eddington accretion of baryon matter, 
the mass growth rate Mbh of a BH is presumed to be pro- 
portional to the black hole mass Mbh, namely 

Mbh = MBH/tsai , (1) 
where tsai is the so-called Salpeter timescale (Salpeter 1964) 

tsai = 7^ — = 3.9 X 10 yr 
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with c, L, iEdd and e being the speed of light, the luminos- 
ity, the Eddington luminosity and the radiative efficiency, 
respectively. For constant e and L/L^dd parameters, the BH 
mass grows exponentially in the form of 



MBH(t) = Mo exp [{t - t„)/tsal] 



(2) 



where Mo is the seed BH mass and to is the initial time 
of accretion. Recent observations suggest L/L-Edd < 1 (e.g., 
Vestergaard 2004; McLurc & Dunlop 2004) and e >0. 1-0.2 
(e.g., Yu & Tremaine 2002; Elvis et al. 2002; Marconi et 
al. 2004). The latest magnetohydrodynamic (MHD) simu- 
lations for disc accretion indicate an e higher than the oft- 
quoted value of ~ 0.1 (e.g., Gammie, Shapiro & McKinney 
2004). A higher value of e tends to increase the Salpeter 
timescale tsai and thus makes the mass growth of a SMBH 
via gas accretion more difficult within a short time (e.g., 
Shapiro 2005). Given an estimated age of ~ 0.9 Gyr for 
the quasar SDSS 1148-1-5251 at 2 = 6.43, it would not be 
easy to assemble a SMBH of mass ~ 3 x 10® Mq from a 
~ 10 — 100 M© seed BH (e.g., the remnant of an imploding 
core of a massive Population HI star; Arnett 1996; Heger 
k. Woosley 2002) even for a sustained accretion of baryon 
matter at the Eddington limit all the time. 

While speculative to various extents, possible theoret- 
ical resolutions to this dilemma of rapid SMBH growth in 

^ In this context, we note in passing the recent detection of a 
gamma-ray burst afterglow with a high redshift z > 6. 



the early Universe include: (1) more massive seed BHs either 
from collapses of supermassive stars (e.g., Shapiro 2004) or 
from accretion of low angular momentum baryon materials 
in the early Universe (e.g., Koushiappas, Bullock & Dekel 
2004); (2) more frequent BH merging processes (e.g., Yoo & 
Miralda-Escude 2004; Shapiro 2005; but see Haiman 2004); 
(3) rapid mass growths via a sustained super-Eddington ac- 
cretion (e.g., Ruszkowski & Begelman 2003; Volonteri & 
Rees 2005). All these proposals with various assumptions 
might produce the required mass > IO^Mq of a SMBH at 
z = 6.43 through a baryon accretion alone. 

Alternatively, a sustained accretion of self-interacting 
dark matter (SIDM) particles (e.g., Spergel & Steinhardt 
2000) onto a seed BH have been modelled to reproduce the 
observed A^bh — o-f, relation (Ostriker 2000; Hennawi & Os- 
triker 2002; cf. MacMUlan & Henriksen 2002 for an alter- 
native approach). As a different application of these ideas 
and as a theoretical contest, here we entertain the possibil- 
ity that a proper combination of SIDM and baryon accretion 
at distinct stages might lead to desired features of forming 
SMBHs in the early Universe. It is natural and sensible to 
imagine that on large scales, SIDM particles and baryons 
are intermingled and mediated by gravitational interactions 
through fluctuations in the early Universe. Based on the the- 
oretical knowledge of accreting baryon matter, we therefore 
advance in this paper a two-phase scenario involving accre- 
tion of both quasi-spherical SIDMs and baryon matter. In 
§2, we describe and elaborate our two-phase accretion model 
scenario in specifics. Summary and discussion are contained 
in §3. 



2 THE TWO-PHASE ACCRETION SCENARIO 

In our two-phase accretion model for SMBH formation, the 
first phase is featured by a sustained, rapid quasi-spherical 
and quasi-steady Bondi accretion of mainly SIDM particles 
(a small fraction of baryon matter mixed therein) onto a 
seed BH created at 2 < 30 presumably by a core implosion 
inside a first-generation massive star of Pop HI until reach- 
ing a BH mass of Mbh ~ 10® Mq during the redshift range 
of z ~ 20 — 15. The second phase of subsequent BH mass 
growth is primarily characterized by a baryon accretion to 
eventually assemble a SMBH of mass Mbh ~ 10® Mq; such 
SMBHs may form either around z ~ 6 for a sustained baryon 
accretion at the Eddington limit or later at lower z for aver- 
age accretion rates below the Eddington limit. For concep- 
tual clarity, we consider these two major phases separately. 
On the theoretical ground, the first phase should gradually 
evolve into a diffusively limited accretion of SIDM parti- 
cles continuing towards the BH. As time goes on, the initial 
mixture of SIDM particles and baryon matter will eventually 
become more or less separated during the accretion process 
in the sense that radiative baryon matter gradually flatten 
to a disc accretion which eventually overwhelms the steady 
accretion of SIDM particles fading into the diffusively lim- 
ited process. 
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2.1 Phase I: a sustained spherical accretion of 

mainly SIDM psirticles onto a seed black hole 

We presume that the initial seed BHs were created by core 
implosions of massive Pop III stars with typical remnant BH 
masses of Mo ~ 10— IOOMq in the redshift range z 30— 10 
[e.g., Wilkinson Microwave Anisotropy Probe (WMAP) ob- 
servations for the excess power in cosmic microwave back- 
ground provide tantalizing evidence for the reionization era; 
see Kogut et al. 2003]; frequent coalescences of such seed 
BHs might possibly happen to produce more massive seed 
BHs around the same epoch or shortly thereafter. When 
such a seed BH happens to immerse in a dark matter (DM) 
halo of low angular momentum, it accretes SIDMs together 
with a small fraction of baryon matter mixed therein. As an 
optimistic approximation, such a SIDM accretion is envis- 
aged as grossly spherically symmetric by an effective trans- 
port of angular momentum outward in the ensemble of 
SIDM particles (Ostriker 2000). We define the specific cross 
section as ctq = Ox/mx for an SIDM particle with a mass 
rux and a cross section Gx; the mean free path is therefore 
A = l/{p(Jo) with p being the mass density of the SIDM in- 
cluding a small mass fraction of baryon matter. For regions 
within the radial range r >\, the SIDM is sufficiently dense 
and may be grossly perceived as a 'fluid' (e.g., Peebles 2000; 
Subramanian 2000; Moore et al. 2000; Hennawi & Ostriker 
2002). 

As a classical reference of estimates, we begin with 
the well-known stationary Bondi (1952) accretion. Given 
a singular isothermal sphere (SIS) mass density profile'^ 
p = C'i / {2-KGr'^) for a SIDM halo, the mass growth with 
time of a BH embedded in a quasi-spherical SIDM halo is 
Mii^{t) = 2Clt/G, where G is the gravitational constant 
and Cs is the SIDM halo "sound speed" (Ostriker 2000). 
Here the "sound speed" is essentially the local velocity dis- 
persion of SIDM particles and is equal to the virial velocity 
for the SIS case. A quantitative comparison of the Bondi 
accretion with the Eddington accretion is given by the ratio 

^Bon t ( Cs Yf AfBH 

MEdd" O.lUokms-iy' WmqJ ' 

where M-boh and MEdd are the Bondi and Eddington mass ac- 
cretion rates, respectively. Apparently, given a seed BH mass 
and a typical halo sound speed (see below) , the Bondi mass 
accretion rate Mboti dominates over the Eddington mass ac- 
cretion rate M-Edd- As SIDM particles do not radiate, this 
super-Eddington accretion will proceed without impedence. 
For sustained isothermal spherical self-similar collapses or 
accretion (e.g., Lou & Shen 2004; Shen & Lou 2004), the 
maximum mass growth rate remains in the same order of 
magnitude as that of the steady Bondi accretion. We em- 
phasize that in the presence of accretion shocks, the central 

While being highly speculative for a SIDM halo, the SIS density 
profile for SIDM and gas was also considered earlier by Ostriker 
(2000) and King (2003). Ostriker discussed consequences of other 
density profiles of r^^ (e.g., Navarro, Frenk, & White 1995) and 
r^'"*/"^ (e.g.. .ling & Suto 2000; Subramanian, Cen, & Ostriker 
2000; Lou & Shen 2004; Bian & Lou 2005; Yu & Lou 2005). 
There are other cusped power-law density profiles available in 
the so-called hypervirial family (e.g., Evans & An 2005). 




t[yr] 



Figure 1. The mass growth history of a black hole in the SIDM 
halo mixed with a fraction of baryon matter. The upper and lower 
panels are for the histories of mass accretion rate dM/dt and of 
black hole mass Mbh growth, respectively. The thin solid line 
before ti represents the Phase I steady Bondi accretion of SIDM 
particles, while the other thin lines are the diffusive limited SIDM 
accretion after ti . The curves of thicker (boldface) linetypes are 
the Phase II baryon accretion at the Eddington rate: here, we 
consider three different cases of the moment that Phase II ax;cre- 
tion begins, that is, before (dotted), at (dash-dotted) and after 
(dashed) ti, respectively. Please note that the dash-dotted line 
and dotted line almost coincide to the right side of ti. The initial 
seed black hole mass due to a Pop III star is 30 Mq at z = 20. 
The input parameters are Cs = 30 km s~^, uq = 1.0 cm'^ g~^ 
and e = 0.15. 



mass accretion rate should be modified (Shen & Lou 2004; 
Bian & Lou 2005). 

This SIDM accretion phase continues until the mean 
free path A becomes comparable to the accretion radius ta 
with a corresponding timescale ti = aoCs/i'^irG), a BH mass 
MBH(ti) = aoCs / {2nG^) and a transitional accretion radius 
rA{ti) = <ToCs/(27rG). For typical parameters, we have the 
following quantitative estimates 

ii ^ 1.1 X lO'' yr (a/30 km s'^) cr„/(l cm^ g-^) , 

Mbh,*! ^ 1.4 X lO'*M0(Cs/3O km s-i)Vo/(l cm^ g"^) , 

rA.ti - 7 pc (C,/30 km s-')Vo/(l cm^g-^) . 

(4) 

For a virialized SIDM halo at high z, the typical halo 
virial velocity or local velocity dispersion of SIDM particles 
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(mimicked as a 'sound speed') may be estimated by 

Cs = 8.2(M/lO^A/0)'/^(l + zf'^ km s"' (5) 

(e.g., Barkana & Loeb 2001). Therefore, an estimate of Cs ~ 
30 km s"^ for a virilized SIDM halo of mass M ~ IO^A/q 
during z ~ 20 — 15 appears plausible. The comoving halo 
number density n{M, z) with mass M at a given z can be 
calculated from the standard hierarchical structure forma- 
tion model. We adopt an input power spectrum computed 
by Eisenstein & Hu (1999). For cosmological parameters, we 
take = 0.3, Oa = 0.7, fib = 0.045, K = 0.7, erg = 0.9 
and n = 1 in our model calculations. 

Estimated constraints on cto from both observations and 
theories are briefly summarized at this point. Wandelt et al. 
(2000) evaluated the constraints on ctq and found a ctq range 
of ~ 0.5 — 6 cm^ g~^. Yoshidaet al. (2000) numerically simu- 
lated the evolution of a galaxy cluster for ctq = 10, 1, and 0.1 
cm^ g~^, and obtained corresponding radial mass profiles. 
Using the high-resolution X-ray data of Chandra satellite 
and the assumption of a hydrostatic equilibrium, Arabadjis 
et al. (2002) derived a mass profile for the galaxy cluster MS 
1358+62 that peaks strongly in the central region. From a 
comparison with the result of Yoshida et al. (2000), they 
concluded that cro < 0.1 cm^ g~^. However, Markevitch et 
al. (2004) pointed out that there are certain difficulties with 
this stringent limit and they provided a less stringent limit 
of cro < 1 eve? g~^. 

In general, ctq may well depend on relative velocity v of 
SIDM particles (e.g., Firmani et al. 2000). Here, we tenta- 
tively adopt a (To in the form of 

/30 km s"^~ 



cm g 



(6) 



where v may be estimated by the isothermal 'sound speed' 
Cs parameter. 

Around and after the time ti, a transition to a slower 
diffusively limited accretion of SIDM particles will gradu- 
ally occur. The accretion rate is then determined by the 
rate at which SIDM particles are scattered into the loss 
cone. Meanwhile, the embedded or intermingled baryon mat- 
ter become more concentrated and flattened through radia- 
tive processes, and their accretion rate towards the central 
black hole becomes more and more important as the dif- 
fusively limited accretion of SIDM particles approaches a 
quasi-steady state. The BH continues to accrete SIDM par- 
ticles further but at a much slower rate than the Eddington 
mass accretion rate of baryons as shown by the mass accre- 
tion rate ratio below from equation (5) of Ostriker (2000) 



^=3.4x10-^ ' 

MEdd 



cro 



0.1 1 cm2 g-i ^30 km s-i 



Cs 



Mbh 
106 A/q 



(7) 

For order-of-magnitude estimates, we may safely ignore the 
SIDM accretion more or less after the formation of a signifi- 
cantly enlarged seed BH around the end of Phase I accretion. 
In this sense and in reference to the very initial seed BH, we 
regard it as a 'secondary seed BH' for the Phase II accretion 
of baryon matter. In Figure^ we have explored three differ- 
ent onset times for the Eddington accretion rate of baryon 
matter for comparison. Within our scenario, it is more sen- 
sible to have the baryon accretion all along with the SIDM 
accretion, roughly corresponding to the heavy dotted curve. 



2.2 Phase II: a disk accretion of baryon matter 

In contrast to baryon matter accretion at the Eddington 
limit as estimated by equation (2), a rapid quasi-spherical 
accretion of SIDM particles during Phase I dominated dur- 
ing the first ~ lO"" yr or so and the BH mass rapidly grows to 
A/bh, ti ~ 10® Mq to become a secondary seed BH for the 
subsequent baryon matter accretion. After this almost in- 
stantaneous Phase-I accretion in reference to the Salpeter 
timescale ts^i, the accumulation of SIDMs proceeds at a 
much slower pace with an inefficient loss cone accretion (Os- 
triker 2000; Hennawi & Ostriker 2002). With favourable en- 
virons or sustained reservoirs of fuels, accretion of baryon 
matter gradually become dominant to increase the BH mass 
by a factor of ~ 10^ within subsequent several Salpeter times 
[see equation (2)] at the Eddington accretion rate. It is this 
subsequent accretion of baryon matter that eventually as- 
sembles most of the SMBH mass, consistent with observa- 
tions that the BH mass densities in nearby galactic bulges 
and in active SMBHs are comparable to the mass density 
accreted during the optically bright/obscured QSO phase 
(e.g., Yu & Tremaine 2002; Fabian 2004; Haehnelt 2004). 
Based on the explorations in Figure 1, it would not matter 
that much as for the exact moment when the Phase II Ed- 
dington baryon accretion sets in. In fact, it may begin at 
any moment (either before or after t\). However, the Phase 
I SIDM accretion can produce a massive enough seed BH 
for the Phase II baryon accretion. 

Figure 1 illustrates a black hole growth history includ- 
ing the two phases. The Phase I SIDM accretion increases 
the black hole mass substantially, and then the Phase II 
baryon accretion enhance the enlarged seed black hole fur- 
ther. The accretion due to SIDM particles after t\ peters 
out rapidly into the diffusively limited regime. 



2.3 Model applications to high-z quasars 

Figure 2 shows the mass range of forming high— 2 SMBHs, 
where we take on three observed SDSS high— 2: quasars with 
reported SMBH masses (Fan et al. 2001, 2003; Willott et al. 
2003; Vestergaard 2004): J114816.64-H525150.3 {z = 6.43, 
ATbh - 3 X 10^ Mq), J103027.10+052455.0 (z = 6.28, 
ATbh - 3.6 X lO'' Mq), and J130608.26+035626.3 {z = 5.99, 
A^BH ~ 2.4 X 10^ Mq) distinguished by dashed, dotted, and 
dash-dotted curves, respectively. Based on our model sce- 
nario, we trace BH masses back to higher z, assuming a 
sustained baryon Eddington accretion and three different ra- 
diative efficiencies e = 0.1, 0.15, 0.2 to compute the required 
minimum BH mass of the phase-I accretion as a function of 
redshift z. We plot the actual final BH mass A^bh, ti of the 
phase-I accretion (light solid curves in Fig. 2) in DM halos 
with 4-0- and 5-o- fluctuations at a given z, using equation 
and Cs given by equation . As an exploration, we also 
calculate A^bh, ti with Cs ^ 100 km s"'^ yet with a smaller 
specific cross section cro — 0.02 cm^ g~^ shown by the heavy 
dashed line in Fig. 2. 

As shown in Fig. 2, our model (i.e., light solid curves) 
can readily account for the three high— z SMBHs for e = 0.1 
and 0.15, without invoking either super-Eddington baryon 
accretion or numerous BH merging processes. For a higher 
e = 0.2, it is unnatural to explain the presence of high- 
z SMBHs with 4-cr fluctuations. For 5— ct fluctuations, the 
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Figure 2. Three sets of dashed, dotted, dash-dotted curves with 
radiative efficiency t = 0.2, 0.15, 0.1 respectively are the required 
minimum BH masses for 10 < z < 30 after the phase-I SIDM ac- 
cretion, with distinct hne types referring to three reported high— z 
quasars (SMBHs). The heavy solid lines show the virialized dark 
matter halo mass Mvir(z) with 4(7 and Scr fluctuations at different 
2 values. The light solid lines show the secondary seed BH mass 
Mbh, ti created by the first rapid accretion phase at various z for 
Cs and (TO according to equations J3 ^'^'^ IHJ- The heavy dashed 
line shows Mbh, ti with a constant sound speed Cs = 100 km s^^ 
and a smaller cross section ctq = 0.02 cm^ g"'^. 

mass requirement can be just met. It is clear that only un- 
der rare circumstances, extremely massive SIDM halos may 
give birth to seed BHs with required Phase-I masses in our 
scenario. 

There are fundamental differences between the two 
main accretion phases, among which the most important one 
is that baryon accretion produces strong detectable feedback 
(e.g., intense radiations and outflows or jets etc) into sur- 
roundings. Reviving some ideas of Silk & Rees (1998), King 
(2003) found a simple yet remarkable association between 
accretion and outflows which plausibly leads to the local 
Man — Ob relation (e.g., Tremaine et al. 2002) at the termi- 
nation of accreting a SMBH. The phase-II accretion in our 
scenario involves processes largely similar to those outlined 
by King (2003) , including Eddington luminosity, baryon ac- 
cretion, intense outflows, and so forth; an Mbh — ot, relation 
(see equation 15 of King 2003) can be established as the 
baryon accretion goes on and the SMBH is continuously as- 
sembled. As described, such an Mbh — Ob relation may not 
emerge before ending the phase-I accretion. To be specific, 
we combine equations (14) and (7) of King (2003) as 

Mbh = ^(^)' ^ 1-5 X 10V^200(«™.M)' Mr, , (8) 

where k is the electron scattering opacity, fg ^ 0.16 is the 
gas fraction, Vm is the mass shell velocity driven by outflows 
and (76200 is the bulge stellar velocity dispersion ab in unit of 
200 km s~^. Potentially, it would also be very important to 
incorporate effects of magnetic field in the baryon accretion 
phase (Yu & Lou 2005; Wang & Lou 2005 in preparation). 
During the accretion phase, Vm < Ob at first and gradu- 
ally Vm approaches ~ ab as the accretion diminishes; mean- 
while, the Mbh — relation appears. During the QSO/AGN 



phase, BH masses may drop below the Mbhu — o"b relation for 
normal galaxies, as indicated by some results of narrow-line 
Seyfert 1 galaxies (e.g., Bian & Zhao 2004; Grupe & Mathur 
2004). Further observations and data analysis of BH masses 
versus velocity dispersions are needed to explore these as- 
pects. 

On the basis of our scenario, we may accommodate the 
high— z and low— 2: quasars in a unified framework: There ex- 
ists a distribution of baryon accretion rates onto secondary 
seed BHs of ~ 10^ Mq produced between 2: ~ 20 — 15 by 
rapidly accreting SIDMs during Phase I. Given specific sit- 
uations, such secondary seed BHs may continue to accrete 
baryon matter with high rates until they run out of fuels 
or the accretion rate becomes considerably lower than the 
Eddington limit; this then gives rise to high— z; SMBHs of 
~ ]^q8-9 Those secondary seed BHs with lower rates of 
accreting baryon matter will evolve into SMBHs much later 
as already discussed above. Given the paucity of high— z 
quasars detected so far, our scenario implies that BH sys- 
tems with high accretion rates should be extremely rare. 



3 SUMMARY AND DISCUSSION 

We summarize the two-phase accretion model of SMBH for- 
mation in the early Universe below and discuss a few obser- 
vational implications. 

Phase I involves a rapid quasi-spherical accretion of 
mainly SIDM particles mixed with a small fraction of baryon 
matter onto an initial seed BH created by a Pop HI star. 
Such a BH grows rapidly to ~ 10® Mq within a fairly short 
timescale of ~ 10^ yr. 

Phase II involves an accretion of primarily baryon mat- 
ter (initially mixed with SIDM particles) at the Eddington 
limit to accumulate most of the BH mass. Since the BH 
mass is sufficiently massive at the beginning of this phase 
II, it takes only several Salpeter e— folding times to grow 
a BH mass of ~ 10^ Mq according to the estimate of an 
exponential growth by equation (2). 

The transition from Phase I to the much slower dif- 
fusion limited accretion of SIDM particles goes on concur- 
rently with the gradual dominance of baryon accretion at 
the Eddington limit in Phase II. 

The Phase-I accretion of SIDMs is crucial to produce 
a sufficiently massive secondary seed BH for further growth 
by accreting baryon matter. In this scenario, the reported 
high- 2 SMBHs of ~ 10^ Mq can be produced, without in- 
voking the hypotheses of either super-Eddington baryon ac- 
cretion or extremely frequent BH merging processes. Based 
on a sample of six quasars with 2 > 5.7 observed by the 
SDSS, Fan et al. (2003) estimated a comoving density of 
such luminous quasars at 2 ~ 6 and found these quasars 
showing a ~ 5(7 peak in the density field. This inference 
may be readily accounted for by our cosmological model re- 
sults shown in Fig. 2. We attempt to combine the models 
of quasi-spherical and quasi-steady SIDM accretion with a 
baryon accretion to give more plausible schemes of accretion 
leading to early formation of SMBHs with masses ~ IO^Mq 
at high redshifts of 2 > 5 — 6. 

We now elaborate on consequences of this two-phase 
scenario. If most SMBHs form by this two-phase accretion 
with fairly rare BH mergers, the currently observed upper 
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mass limit for a SMBH < 10^ Mq constrains the combina- 
tion of the specific SIDM cross section ao and the effective 
baryon accretion time t2 (cf. equation 4), such that 

, ? -1 f ,nn^° -0 ' e^P(*2 AsaO < 60 , (9) 
f cm g V 100 Km s J 

where the Eddington luminosity is adopted.^ For SIDM ha- 
los formed at low z, they may become massive enough to 
make the virialized velocity dispersion or 'sound speed' of 
the order of Cs ~ 100 km s~^. So the baryon accretion 
time t2 cannot exceed ~ Stsai for low— z QSOs. For typi- 
cal parameters of e = 0.1 and L — L^dd, it would require 
a t2 < 2 X 10* yr, consistent with the observational QSO 
lifetime of tq ~ 10^ — 10* yr estimated for low— z QSO 
populations (e.g., Yu & Tremaine 2002; McLure & Dunlop 

2004) . 

If we take on reference values of equation @ for the 
particular z = 6.43 SMBH with ctq = 1 cm^ g~^ and 
Cs ~ 30kms~^, it would have accreted baryon matter 
over ~ Stsai (~ 3 X 10* yr) for e = 0.1 and L = L^dd, 
in the absence of accretions at super-Eddington rates and 
of BH merging processes. Such a Phase II accretion time 
is much longer than the low— 2; QSO lifetime, but accounts 
for only ~ 1/3 of ages for these quasars. We speculate that 
high— quasars may have longer effective accretion times 
than low— 2 quasars as the former might have more abun- 
dant fuel supplies. The plausible physical reasons include: 
(1) gas materials in galaxies have been much consumed by 
star formation activities in the low— 2 quasars; (2) interac- 
tions among galaxies might have been more frequent in the 
early Universe, that may trigger high accretion activities. In- 
sufficient or interrupted baryon accretion would very likely 
leave behind less massive SMBHs with masses well below 
the Mbh — o"6 relation curve. 

We note that the Phase I accretion will not affect ob- 
servations at low— 2 galaxies. First, the Phase I SIDM Bondi 
accretion is more dominant than the Eddington accretion of 
baryons only for very early {z > 10) SIDM halos, as such ha- 
los tend to have a high 'sound speed' Cs (cf. equations 3 — 5). 
Secondly, the phase-I accretion might be severely limited by 
the inner density profile of the SIDM halo. We have used a 
SIS mass density profile p oc r~^. For a self-similar accretion 
at a given time t, one would have p oc r~^^'^ (e.g., Jing & 
Suto 2000; Lou & Shen 2004; Bian & Lou 2005; Yu & Lou 

2005) . Numerical N-body simulations (e.g., Navarro et al. 
1997) indicate an inner p cx r^^ , i.e. the NFW profile (see 
discussions of Ostriker 2000) . For the more inclusive family 
of the hypervirial models, we have p cx r"'^"^' with the in- 
dex p falling in the range of < p < 2 (Evans & An 2005) . 
A shallower profile may lead to a less efficient accretion of 
SIDMs (e.g., Hennawi & Ostriker 2002). Observations find 
that the inner density profiles (e.g., a p oc r~°'^) of dark 
matter halo of nearby galaxies are shallower than both the 
SIS and NFW density profiles, and therefore the resulting 
secondary seed BH mass A^bh, ti may be expected to be 
smaller. 



Should we take into account of BH merging processes by a sig- 
nificant mass amplification factor of ~ 10* (e.g., Yoo & Miralda- 
Escude 2004) , then this upper limit would be reduced by a factor 
of ~ 10* accordingly. 



Specific to our Milky Way galaxy, the central black hole 
has been inferred to possess a mass of ~ 4 X 10*^ Mq by steUar 
dynamic diagnostics. For such a less massive SMBH, there 
are many possible ways to assemble such a BH and it is not 
necessary to invoke our two-phase scenario. 

Finally, we note possible observational signatures and 
consequences of the two-phase scenario envisaged in this pa- 
per. If some large regions in the early Universe were some- 
how distributed with considerable less baryon matter (e.g., 
the gravitational potential well of the halo may not be deep 
enough to bind the high speed baryons due to supernova 
or hypernova feedbacks), then we may find some DM halos 
with BHs at their centres but without forming host galaxies. 
In other words, these BHs grow entirely by Phase I SIDM 
quasi-spherical Bondi accretions ended with diffusively lim- 
ited phase, and the baryon accretion has never occurred in 
a significant manner. Such unsual systems of DM halo-BHs 
might be revealed directly by chance through gravitational 
lensing effects. 

Very recently, Magain et al. (2005) reported the discov- 
ery of a bright quasar HE0450-2958 yet without a massive 
host galaxy. The black hole of HE0450-2958 may well be em- 
bedded within a dark matter halo ('dark galaxy'), constitut- 
ing a DM halo-BH system. Ambient interactions with such a 
massive object could more readily explain the ring-like star- 
burst in the neighbouring galaxy as well as the capture of 
gas materials, leading to the eventual onset of quasar ac- 
tivities we observe. Our model can readily account for such 
kind of phenomena. The phase I SIDM accretion in the dark 
matter halo produced a seed black hole with a mass range 
of ~ 10^ — 10^ Mq, but the phase II baryon accretion has 
never occurred due to the absence of a host baryon galaxy. 
As such a massive seed black hole travelled across a neigh- 
bouring disc galaxy, it began to induce baryon accretion 
activities violently to give rise to a bright quasar. We will 
further develop models on the basis of such a scenario in a 
separate paper. 

Another interesting piece of observational evidence 
comes from the galaxy NGC 4395 (Peterson et al. 2005). 
A very 'smaU' low-luminosity SMBH (~ 3 x 10^ Mq) re- 
sides in the bulgeless galaxy NGC 4395, implying that the 
stellar bulge is not a necessary prerequisite for a black hole 
in the nucleus of a galaxy. Our model can provide a plausible 
explanation in the sense that a relatively small SMBH is the 
product of the phase I SIDM accretion, while the phase II 
baryon accretion at the Eddington rate is almost absent in a 
bulgeless galaxy (of course, a very low rate baryon accretion 
might take place). 
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